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PhI Abstract 

rr^' The method for simulation of medium-induced rescattering and energy loss of hard 

r~| , partons in ultrarelativistic heavy ion collisions is developed. The model is realized as fast 
Monte-Carlo tool implemented to modify standard PYTHIA jet event. 
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1 Introduction 

The experimental investigation of ultra-relativistic nuclear collisions offers a unique possibility 
of studying the properties of strongly interacting matter at high energy density. In that regime, 
hadronic matter is expected to become deconfined, and a gas of asymptotically free quarks and 
gluons is formed, the so-called quark-gluon plasma (QGP), in which the colour interactions 
between partons are screened owing to collective effects [1]. One of the important tools to study 
QGP properties in heavy ion collisions is a QCD jet production. Medium-induced energy loss 
of energetic partons, the so-called jet quenching, has been proposed to be very different in cold 
nuclear matter and in QGP, resulting in many challenging observable phenomena [2]. Recent 
RHIC data on suppression of inclusive high-p^ charge and neutral hadron production from 
STAR [3], PHENIX [4], PHOBOS [5] and BRAHMS [6] are in agreement with the jet quenching 
hypothesis [7]. However direct event-by-event reconstruction of jets and their characteristics is 
not available in RHIC experiments at the moment, while the assumption that integrated yield 
of all high-pT particles originates only from jet fragmentation is not obvious. 

At LHC a new regime of heavy ion physics will be reached at ^snn = 5.5 TeV where 
hard and semi-hard QCD multi-particle production can dominate over underlying soft events. 
The initial gluon densities in Pb— Pb reactions at LHC are expected to be significantly higher 
than at RHIC, implying stronger partonic energy loss which can be observable in various new 
channels [8]. Thus in order to test a sensitivity of accessible at LHC observables to jet quenching, 
and to study corresponding experimental capabilities of real detectors, the development of fast 
Monte-Carlo tools is necessary. 

2 Physics frameworks of the model 

The detailed description of physics frameworks of the developed model can be found in a number 
of our previous papers [9, 10, 11, 12, 13]. The approach relies on an accumulative energy losses, 
when gluon radiation is associated with each scattering in expanding medium together including 
the interference effect by the modified radiation spectrum dE/dl as a function of decreasing 
temperature T. The basic kinetic integral equation for the energy loss Ai5 as a function of 
initial energy E and path length L has the form 

^ ' 

where / is the current transverse coordinate of a parton, dP/dl is the scattering probability 
density, dE/dl is the energy loss per unit length, A = l/{ap) is in-medium mean free path, 
p oc T^ is the medium density at the temperature T, a is the integral cross section of parton 
interaction in the medium. 

The collisional energy loss due to elastic scattering with high-momentum transfer have been 
originally estimated by Bjorken in [14], and recalculated later in [15] taking also into account the 



loss with low-momentum transfer dominated by the interactions with plasma collective modes. 
Since latter process contributes to the total collisional loss without the large factor ~ In {E/fiu) 
{fiD is the Debye screening mass) in comparison with high-momentum scattering and it can 
be effectively "absorbed" by the redefinition of minimal momentum transfer t ~ /i|, under the 
numerical estimates, we used the collisional part with high-momentum transfer only [10], 
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and the dominant contribution to the differential cross section 
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for scattering of a parton with energy E off the "thermal" partons with energy (or effective 
mass) rriQ ~ 3T ^ E. Here C = 9/4,1,4/9 for gg, gq and qq scatterings respectively, ag 
is the QCD running coupling constant for Nf active quark flavors, and Aqcd is the QCD 
scale parameter which is of the order of the critical temperature, Aqcd — T^ — 200 MeV. 
The integrated cross section a is regularized by the Debye screening mass squared fJ^j^iT) ~ 
47rQ;sT^(l + A^j/6). The maximum momentum transfer tinax = [s — {iTip+mQy][s — {mp~mQy]/s 
where s = 2mQE + ml + rrip, rrip is the hard parton mass. 

There are several calculations of the inclusive energy distribution of medium-induced gluon 
radiation from Feyman multiple scattering diagrams. The relation between these approaches 
and their main parameters have been discussed in details in the recent writeup of the working 
group "Jet Physics" for the CERN Yellow Report [8] . We restrict to ourself here by using BDMS 
formalism [16]. In the BDMS framework the strength of multiple scattering is characterized by 
the transport coefficient q = f^j^/Xg {\g is the gluon mean free path), which is related to the 
elastic scattering cross section a (3). In our simulations this strength in fact is regulated mainly 
by the initial QGP temperature Tq. Then the energy spectrum of coherent medium-induced 
gluon radiation and the corresponding dominated part of radiative energy loss of massless 
parton has the form [16]: 
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where ri = L/{2\g), y = uj/E is the fraction of the hard parton energy carried by the radiated 
gluon, and Cr = 4/3 is the quark color factor. A similar expression for the gluon jet can be 
obtained by substituting Cr = 3 and a proper change of the factor in the square bracket in (4), 
see ref. [16]. The integral (4) is carried out over all energies from ij„im = E^pm = l^'o^g, the 
minimal radiated gluon energy in the coherent LPM regime, up to initial jet energy E. 



The generalization of the formula for heavy quark of mass m^ was done by using "dead-cone" 
approximation [17]: 

—\ 1 dE ] - ( A.y^^ f!!h.Y^^ 

but note there are exist more recent developments on heavy quark energy loss in the litera- 
ture [18, 19]. 

The medium was treated as a boost-invariant longitudinally expanding quark-gluon fluid, 
and partons as being produced on a hyper-surface of equal proper times r [20]. In order to 
simplify numerical calculations in original version of the model we omit the transverse expansion 
and viscosity of the fluid using the well-known scaling Bjorken's solution [20] for temperature 
and density of QGP at T > T^ ~ 200 MeV: 

eiry/' = e,4'\ Tiry^' = T,rl'\ p{t)t = p.r,. (7) 

For certainty we used the initial conditions for the gluon-dominated plasma formation expected 
for central Pb— Pb collisions at LHC [21]: 

To ~ 0.1 fm/c. To ~ 1 GeV, Pg ^ l.95T^ . 

Then for non-central collisions and for other beam atomic numbers we suggest the propor- 
tionality of the initial energy density Eq to the ratio of nuclear overlap function and effective 
transverse area of nuclear overlapping [10]. 

Note, however, that using other initial parameters and scenarious of QGP space-time evolu- 
tion for Monte-Carlo realization of the model is possible (by changing some internal parameters 
of the routine). In fact, the influence of the transverse flow, as well as of the mixed phase at 
T = Tc, on the intensity of jet rescattering (which is a strongly increasing function of T) has 
been found to be inessential for high initial temperatures Tq ^ Tc. On the contrary, the pres- 
ence of QGP viscosity slows down the cooling rate, which leads to a jet parton spending more 
time in the hottest regions of the medium. As a result the rescattering intensity goes up, i.e., 
in fact an effective temperature of the medium gets lifted as compared with the perfect QGP 
case. We also do not take into account here the probability of jet rescattering in nuclear matter, 
because the intensity of this process and corresponding contribution to total energy loss are 
not signiflcant due to much smaller energy density in a "cold" nuclei. 

Another important set of the model is the angular spectrum of in-medium gluon radiation. 
Since the full treatment of angular spectrum of emitted gluons is rather sophisticated and 
model-dependent [9, 16, 22, 23, 24], simple parameterizations of gluon angular distribution 
over emission angle 6 was used: 

ocsm^exp — , (8) 
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where d^ ~ S'' is the typical angle of coherent gluon radiation estimated in work [9]. Other 
parameterizations are also possible. 



3 Monte-Carlo simulation procedure 

The model has been reahzed as fast Monte-Carlo event generator, and corresponding Fortran 
routine PYQUEN.F is available by the web [25]. The following input parameters should be 
specified by user to fix QGP properties: beam and target nucleus atomic number and type 
of event centrality selection (options "fixed impact parameter" or "minimum bias events" are 
foreseen). Since the routine is implemented as a modification of standard PYTHIA6.2 jet 
event [26], the main user program should be compiled with this version of PYTHIA. 
The following event-by-event Monte-Carlo simulation procedure is applied. 

• Generation of initial parton spectra with PYTHIA (fragmentation off). 

• Generation of jet production vertex at impact parameter b according to the distribution 
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where ^i 2(&, t, ■?/') are the distances between the nucleus centers and the jet production 
vertex \/(r cos ■?/',r sin ■?/'); rmax{b,ip) < Ra is the maximum possible transverse distance r 
from the nuclear collision axis to the V; Ra is the radius of the nucleus A; T4(ri 2) is the 
nuclear thickness function (see ref. [10] for detailed nuclear geometry explanations). 

Generation of scattering cross section da/dt (3). 

Generation of transverse distance between scatterings, /j = (tj+i — Ti)E /pt'- 
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Reducing parton energy by coUisional (2) and radiative (4), (6) loss per scattering i: 

Calculation of parton transverse momentum kick due to elastic scattering i: 

Akl = {E- -^f -{p- --^ - ^f - ml 
2moi p 2moi 2p ^ 

Halting parton rescattering if 1) parton escapes from dense zone, or 2) QGP cools down 
to Tc = 200 MeV, or 3) parton loses so much energy that its Pt{t) drops below 2T(r). 

In the end of each event adding new (in-medium emitted) gluons into PYTHIA parton 
list and rearrangements of partons to update string formation are performed. 

Formation of final hadrons by PYTHIA (fragmentation on). 



As the example, let us compare liadron pT^-spectra obtained with jet quenching (PYTHIA+ 
PYQUEN) and without one (PYTHIA only). Fig.l shows this spectrum for ^Spp = 5.5 TeV 
and QGP parameters selected for central Pb— Pb collisions, the events being triggered by having 
at least one jet with Et > 100 GeV in the final state. As it could be expected, energy loss 
of hard partons results in the high-pT suppression of the spectrum, while the intensive gluon 
emission of soft and semi-hard gluons governs the low-p^ enhancement. 
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Figure 1: The hadron pT-spectrum at WSpp = 5.5 TeV obtained with jet quenching 
(PYTHIA+PYQUEN, dashed histogram) and without one (PYTHIA, sohd histogram). QGP pa- 
rameters were selected for central Pb— Pb collisions. Events were triggered having at least one jet 
with Et > 100 GeV in the final state. 



4 Conclusions 

The method to simulate rescattering and energy loss of hard partons in ultrarelativistic heavy 
ion collisions has been developed. The model is realized as fast Monte-Carlo tool implemented 
to modify standard PYTHIA jet event. Corresponding Fortran routine is available by the web. 
To conclude, let us discuss the physics validity of the model application. 

• Internal parameters of the routine for initial conditions and space-time evolution of quark- 
gluon plasma were selected as an estimation for LHC heavy ion beam energies. The result 
for other beam energy ranges, obtained without additional internal parameters adjusting, 
is not expected to be reasonable. 

• Hydro-type description of expanding quark-gluon plasma used by the model can be ap- 
plicable for central and semi-central collisions. The result obtained for very peripheral 
collisions {b ~ 2Rj^) can be not adequate. 

• Physics model for medium-induced gluon radiation is valid for relatively high transverse 
momenta of jet partons (^ 1 GeV/c). Thus setting reasonably high value of minimum 
Pt in initial hard parton sub-processes in PYTHIA is preferable. 
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